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Abstract. We study the Cauchy problem associated with the system of two conservation laws 
arising in isothermal gas dynamics, in which the pressure and the density are related by the 
7-law equation p(p) ~ p 1 with 7 = 1. Our results complete those obtained earlier for 7 > 1. 
We prove the global existence and compactness of entropy solutions generated by the vanishing 
viscosity method. The proof relies on compensated compactness arguments and symmetry 
group analysis. Interestingly, we make use here of the fact that the isothermal gas dynamics 
system is invariant modulo a linear scaling of the density. This property enables us to reduce 
our problem to that with a small initial density. 

One symmetry group associated with the linear hyperbolic equations describing all entropies 
of the Euler equations gives rise to a fundamental solution with initial data imposed to the line 
p = 1. This is in contrast to the common approach (when 7 > 1) which prescribes initial data 
on the vacuum line p = 0. The entropies we construct here are weak entropies, i.e. they vanish 
when the density vanishes. Another feature of our proof lies in the reduction theorem which 
makes use of the family of weak entropies to show that a Young measure must reduce to a Dirac 
mass. This step is based on new convergence results for regularized products of measures and 
functions of bounded variation. 



1. Introduction 

We consider the Euler equations for compressible fluids 

(1.1) d t p + d x ( P u) = 0, 

(1.2) d t {pu)+d x (pu 2 +p(p))=0, 

where p > denotes the density, u the velocity, and p{p) > the pressure. We assume that the 
fluid is governed by the isothermal equation of state 

(1.3) P(p) = k 2 p, 

where k > is a constant. Observe that the scaling u — > ku, t — > t/k allows one to reduce the 
system (1.1)— (1.3) to the same system with k = 1. 

The existence of weak solutions (containing jump discontinuities) for the Cauchy problem 
associated with (1.1)— (1.3) was first established by Nishida [23] (in the Lagrangian formulation). 
The solutions obtained by Nishida have bounded variation and remain bounded away from the 
vacuum. For background on the BV theory we refer to [6l [16] . 

By contrast, we are interested here in solutions in a much weaker functional class and in solu- 
tions possibly reaching the vacuum p — 0. Near the vacuum, the system (1.1)— (1.3) is degenerate 
and, in particular, the velocity u can not be defined uniquely. Indeed, the present paper is devoted 
to developing the existence theory in a framework covering solutions satisfying 

pei°°(n), P \u\ <c(p + p|io g/0 |), n = Rx(o,T), 

with a constant C > depending solely on initial data. The time interval (0, T) is arbitrary. Our 
proof extends DiPerna's pioneering work [10] concerned with the pressure law p(p) ~ p 1 . 
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2. Main result 

Introducing the momentum variable to := pu, one can reformulate the Cauchy problem asso- 
ciated with (1.1)— (1.3) as follows: 

d t p + d x m = 0, 

,2 



(2.1) 



777 

d t m + dJ h p) = 0, 

P 



with initial condition 

(2.2) p\t=o = Po, m\ t =o = m Q := p u . 

where Po,uq are prescribed. Let us first recall the following terminology. A pair of (smooth) 
functions 77 = n(m,p), q = q(m, p) is called an entropy pair if, for any smooth solution (to, p) of 
(|2.1[) . one also has 

d t n(m, p) + d x q(m, p) = 0. 
More precisely, we consider entropies 77,(7 G C 2 (f2) n C 1 (f2) in any domain of the form 

tt := {0 < p < p*, |to| < c*p(l + I In p I ) } , c* > 0, p* > 0. 
It is easily checked that 77, (7 must solve the equations 

2 
2" 



vn in 

(2.3) g m = 2 — 77 m + ?7 P , q p = r\ rn -77™, 

P P 



which implies that 

(2.4) Vpp o 'Huu ttT] uu . 

p z p z 

A pair (77, q) is said to be a weak entropy if 77(0, 0) = (7(0, 0) = 0. It is said to be convex if in 
addition, 77 is convex with respect to the conservative variables (p, m). 
Given an initial data too, po <G L°°(M) obeying the inequalities 

(2.5) po(x)>0, \m (x)\<c po(x)(l + \\np (x)\), xeR 

for some constant Co > 0, an entropy solution to the Cauchy problem (|2.1|) - (|2.2p on the time 
interval (0, T) is, by definition, a pair of functions (to, p) 6 i°°(II) satisfying the inequalities 

(2.6) p(x,t)>0, \m{x,t)\ <cp(x,t)(l + \\np{x,t)\), (x,t)eU 
for some positive constant c, together with the inequality 

(2.7) J J (77(777, p) d t <f + q{m, p) d x ipj dxdt + J ri(m , p ) <p(-, 0) dx > 

for every convex, weak entropy pair (rj,q) and every non-negative function ip € D(K x [0,T)) 
(smooth functions with compact support). 

The main results established in the present paper are summarized in Theorems 2.1-2.3 below. 

Theorem 2.1. (Cauchy problem in momentum-density variables.) Given an arbitrary time in- 
terval (0,T) and an initial data (too,Po) £ L°°(R) satisfying the condition (|2.5j) . there exists an 
entropy solution (to, p) of the Cauchy problem (|2.1[) - (I2.2[) satisfying the inequalities (|2.6[) . with a 
constant c depending on cq only. 

To prove this theorem it will be convenient to introduce the Riemann invariants W and Z by 

W := pe u , Z := pe~ u , 

or equivalently 

p = h{W, Z) (WZ) 1 / 2 , pu = f 2 (W, Z) := (WZ) 1 ' 2 \n(W/Z)^ 2 . 
One can then reformulate the Cauchy problem (|2.ip - (|2.2[) in terms of W, Z , as follows 
d t f 1 (W,Z)+d x f 2 (W,Z) = 0, 

(2 8) 

dth(W,Z)+d x (f 3 {W,Z) + f 1 (W,Z)) = 0, h^iWZf^iXniW/Z) 1 ' 2 ) 2 , 
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(2.9) W\ t=0 = W := Poe 110 , Z\ t=0 = Z := Po e""°. 

a pair of non-negative functions W, Z € L°°(II) if then called an entropy solution to the problem 
if 



few, Z)) ftp + q(W, Z) d x <p) dxdt + / rj{W , Z ) <p{; 0)dx>0 

for any non-negative function <p G D(R x [0, T)), where 

r?(W, Z) := 77(/ 2 (W, Z), A(W, Z)), Z) := g(/ 2 (W, Z), f x {W, Z)), 

and (?7, g) is any convex, weak entropy pair in the sense introduced above. 
Theorem 2.1 above will be obtained as a corollary of the following result. 

Theorem 2.2. (Cauchy problem in Riemann invariant variables.) Given non-negative functions 
Wq,Zq £ L°°(M), the Cauchy problem (|2.8p - (|2.9p has an entropy solution on any time interval 
(0,T). 

It is checked immediately that, if (W, Z) is an entropy solution given by Theorem 2.2, then 
the functions m := f2(W,Z) and p :— f\{W,Z) determine an entropy solution of the problem 

One more consequence of Theorem 2.2 concerns the original problem (|l.l|l - (|1.3|) in the density- 
velocity variables. Defining the density and velocity from the Riemann variables by 

u := ln(WyZ) 1 / 2 , p-^iWZ) 1 / 2 , 

we deduce also the following result from Theorem 2.2. 

Theorem 2.3. (Cauchy problem in velocity-density variables.) Let (0,T) be a time interval. 
Given any measurable functions uq and po satisfying the conditions 

< po e £°°(R), |uo(s)| <co(l + |lnpo(aO|), 

for some positive constant cq, there exist measurable functions u — u(x,t) and p = p{x,t) such 
that 

o<pei°°(n), \u(x,t)\ < c (i + \\n P (x,t)\), (i,t)en 

(where c > is a constant depending on cq) and (u, p) is an entropy solution of the problem 
(11 . 1[) - (j 1 . 3[) in the sense that the entropy inequality 



(r](p,pu)d t ip + q(p,pu)d x ip) dxdt+ / n(p , p u Q ) ip(; 0) dx > 
holds for any convex, weak entropy pair (jj,q) and any function (p as in Theorem 2.1. 



The novel features of our proof of the above results are : 

• the use of symmetry and scaling properties of both the isothermal Euler equations and 
the entropy-wave equation, 

• an analysis of new nonconservative products of functions with bounded variation by mea- 
sures. 

We rely on two classical ingredients. The first tool is the compensated compactness method 
introduced by Tartar in [32J, [33]. (See also Murat [22].) This method allows to show that a 
weakly convergent sequence (of approximate solutions given by the viscosity method) is actually 
strongly convergent: such a result is achieved by a "reduction lemma" (to point mass measures) 
for Young measures representing the limiting behavior of the sequence. Tartar method was applied 
to systems of conservation laws by DiPerna [9[ [10] . For a completely different approach to the 
vanishing viscosity method, we refer to Bianchini and Bressan [2]. Still another perspective is 
introduced in LeFloch [17] , 

The second main tool is the symmetry group analysis of differential equations which goes 
back to Lie's classical works. The first symmetry property we use concerns the system (1.1) 
(1.3) itself: we observe that it is invariant with respect to the scaling p — > Xp (A being an 
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arbitrary parameter). This property allows us to assume that the density is sufficiently small 
when performing the reduction of the Young measures. 

To generate the class of weak entropies, we calculate all the Lie groups associated with the 
entropy equation f|2.4|) for the function r\. By using one of them we construct the fundamental 
solution with initial data prescribed on the line p = 1. This is in contrast with the standard 
approach which prescribes initial data on the vacuum line p = 0. 

The need of a large family of weak entropies for the Young measure reduction was demonstrated 
by DiPerna for the isentropic gas dynamics equations with the pressure law p = p 1 , 7 > 1. When 
7 = 2n+i ' w ith n being integer, DiPerna used weak entropies which are progressive waves given 
by Lax. The method of Tartar and DiPerna was then extended by Serre [29] to strictly hyperbolic 
systems of two conservation laws, by Chen, et al. [3l [8] to fluid equations with 7 G (1,5/3] and 
by Lions, Perthame, Souganidis, and Tadmor [T8l |T9] to the full range 7 > 1. The theory was 
extended to real fluid equations by Chen and LeFloch [H [T5l [5] . We also mention the important 
work by Perthame and Tzavaras on the kinetic formulation for systems of two conservation laws; 
see [26l [27] . The success of these works relies on a detailed analysis of the fundamental solution 
of the entropy wave equation (|2.4[) , which is a degenerate, linear wave equation. 

When 7=1 the analysis developed in [4, 5 for the construction of entropies does not work 
because the equation ()2.4j) degenerates at a higher degree and the Cauchy problem at the line 
p = becomes highly singular. One novelty of the present paper is to rely on symmetry group 
argument to identify the entropy kernel. 

For the convenience of the reader we summarize now the main steps of the proof of Theorems 
2.1-2.3. 

Step 1. We rely on the vanishing viscosity method and first construct a sequence of approximate 
solutions (u e ,p c ), e J, 0, defined on the strip II, and such that 

2e r <p £ <p 2 <l 

for some r > 1. The constant pi can be chosen to be arbitrary small by introducing a rescaled, 
initial density Xp$. We will thus establish first Theorems 2.1 to 2.3 in the case when the initial 
density is small. Then we will treat the general case by observing that the system (1.1)— (1.3) is 
invariant via the symmetry (u, p) — > (it, Xp). More precisely, given an entropy solution (u, p) of 
the problem (1.1)-(1.3) with initial data (uo,po), the pair (v! , p') := (u, Xp) is also an entropy 
solution with the initial data (u' Q , p' Q ) — (ito, Xpo). 

Step 2. Next, we prove that there is a sequence e | such that 

(2.10) W e := p'e u ' W, Z e := p e e' u ' Z weakly * in L^ C {U) 

and there exist Young measures i/ x> t, associated with the sequence e J, and defined on the 
(W, Z)— plane for each point (x, t) £ II, such that 

limF(W e (x,t),Z*(x,t)) = jjF{a,l3)dv Xtt = {v Xtt ,F)=:{F) in D'(M 2 ) 

for any F(a,/3) G C; oc (R 2 ). The crucial point in the compensated compactness argument is to 
prove that v is a point mass measure. In that case the convergence in (I2.10P becomes strong in 
any L[ oc (n), 1 < r < 00. 

Step 3. Given two entropy pairs (rfi, qi), obeying the conditions of Theorem 2.1, we check that 
Tartar's commutation relations 

(2-11) {vx,t,ViQ2-il2Qi) = (vx,t,r)i) { v x,t,q2) - {vx,t,m) {vx.uQi) 

hold. Here, we apply the so-called div-curl lemma of Murat [22] and Tartar [321 S3- The objective 
is to prove that the measure v is a point mass measure by using a "sufficiently large" class of 
entropy pairs in (|2.1ip . 
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Step 4. To produce a large family of entropy pairs, we have to construct a fundamental solution 
x(R,u — s) (where R := lnp) of the entropy equation (|2.4[) . To this end, we rely on symmetry 
group arguments for the equation (|2.4[) . We find that it has an invariant solution 



r,(u,p) = ^pf(u 2 -In 2 p), where £ /"(£) + /'(£) + /(£) = 0. 
Then we define 

X = e fi / 2 /(|w-s| 2 -i? 2 )l|n-s|<|fl|. 
The function /(£) can be represented by a Bessel function of zero index. 

Step 5. Then we search for the entropy pairs in the form 
(2.12) r, = J X (R,u-sMs)ds, q = J a(R,u, s) ds, 

where iji £ L (R) is arbitrary and we describe properties of the kernels Xi a - I n particular, we find 
that a — u\(R,u — s) + h(R,u — s), where the function h is given by an explicit formula. We 
also will show that 

Px ■■= d s x = e R/2 {8s=u-\r\ ~ 8 s =u+\R\) +G x (R,u- s) 1|„_ S |<| R |, 

Ph := d s h = e R/2 U s=u _|ij| + 5 s=u+ \ R ^j + G h (R,u- s) l| u _ a |<|jj|, 
where G x (R,v) and G h (R,v) are bounded, continuous functions. 

Step 6. Finally, we plug the entropy pairs (|2.12[) in Tartar's commutation relations, but in the 
form derived by Chen and LeFloch [3] . We arrive after cancellation of ip at the following equality 
in D'(K) 3 

{Xl P 2h% ~ h 1 P 2 X2)(P3X3) + (/U-P3X3 - XlP?> hl i){P2X2) = -(P3h 3 P2X2 ~ ^3X3^2^2) (Xl) , 

where the notations gi := g(R, U, s,) and Pigi := d Si g(R, u, s,) are used. Then we test this equality 
with the function 

where ip E D(U.) and <pj are molifiers such that 

Y := I / (^2(52) ^3(53) - ^3(^2) ^2 (S3)) ds 2 ds3 ^ 0. 



— oo t/ — oo 



This identity involves products of measures by functions of bounded variation. Such products 
were earlier discussed by Dal Maso, LeFloch, and Murat [7]. 
By letting S go to zero we obtain the equalities 



(2.13) Y D(p)p Jp' dv(W,Z')dv(W,Z)=0, 

JJw,z J J {w<W}n{z'<i/w j 

(2.14) Y ff D(p)p[[ Jp y dv(W',Z')dv(W,Z)=0, 

J Jw,z J J { W'<l/Z \r\\Z' <z\ 



where 



{w<i/z}n{z'<z} 
p = {WZfl 2 , D(p) = V/5(-5 + ylni), p' = (W'Z'Y /2 , 



and the measure du(W , Z') is a copy of dv on the (W , Z')-plane. At this point we choose 
the constant p 2 (see Step 1) small enough to ensure the inequality D(p) > J~p/2- Hence, it 
follows from (j2. 13[) . ()2. 14() that dv x ,t = aSp + p x .t and a (1 — a) = 0, where P(x,t) is a point 
on the (W,Z)— plane and the support of the measure supp p x .t lies in the set{,o = 0}. This 
representation formula for the measure v x , t enables us to justify the passage to the limit as e J. 0. 
We summarize Step 6 in the following key result. 
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Theorem 2.4. Let (m n , p n ) be a bounded in L°°(TV) sequence of entropy solutions of the problem 
(|2.ip and such that 

0<pn, \m n \ < cp n (l + \ln p n \) 
uniformly in n. Then, passing to a subsequence if necessary, (m n , p n ) converges almost everywhere 
in II to an entropy solution (m, p) of (|2 . 1[> . 



3. Vanishing viscosity method 
Given parameters e, t\ > we consider the Cauchy problem 

(3.1) p t + {pu) x = ep xx + 2eiu x , 

(3.2) (pu) t + {pu 2 ) x + p x = e(pu) xx + ei(u 2 ) x + 2ei(lnp) x , 
with initial condition 

(3.3) p\ t= Q = Po + 2ei, u| t=0 = u%. 

In this section we establish the existence of smooth solutions to this problem. Later in this section 
we will assume that t\ = e r for some r > 1. The positivity of the density will be obtained by the 
following argument. 

Lemma 3.1. (Positivity for convection-diffusion equations.) If v — v(x,t) is a smooth bounded 
solution of the Cauchy problem 

(3.4) v t + (uv) x = ev xx , v\ t=0 =v (x), 
where u = u(x,t) £ L°°(n) and uq € L°°(R) 7 then v > provided vq > 0. 

Proof. Given R > 0, let ip : R + — > R be a non- increasing function of class C 2 such that ip( x ) = 1 
for x £ [0,R], ip(x) — e~ x for x > 2R, and ip(x) is a non-negative polynomial for R < x < 2R. 
Denote *(x) = ip(\x\) for i£l Clearly, 

(3.5) |*' (1)1 < \*"(x)\<^(x) 
for some constant c\ > 0. The map 



^ W ~ \ 0, « > 0, 

is a regularization of the mapping v i— y := max{— w, 0}. 

Using (|3.4[) and (|3.5p we can compute the i-derivative of the integral J ^U^v) dx: 

d 2 U u , . _ . f 8U a 

di 

d 2 U u , /• ,dU u 



( 3 - 6 ) = / 1 ^ i ^(e* a; + u*)rfa;+ / v-^(e^ xx + u V x ) dx 



< I ^£v\v x \V(e Cl /R+\u\)dx+ v\^\y{e Cl /R 2 + \u\ Cl /R)dx. 



Observe that 



evi - v\v x \{ec x /R + \u\) = e(\v x \ - + ^)) 2 - v 2 (^ + ^) 2 , 

^ d 2 U, < p 2 v 2 



dv 2 ~ {v 2 +p 2 ) 3 / 2 ' 

and 



dU, 



v— ► V- as it — * 0. 

ov 

We integrate ()3.6|) with respect to t and let /i tend to zero, by taking into account that U^vq) = 0: 

t 



J^v. dx< J J 'i& v _(y* + \ASL)dxdT. 
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By Gronwall's lemma, J dx — 0. We thus conclude that v > 0. □ 

As a consequence of Lemma 3.1, we deduce that any bounded solution (u, p) of the problem 
(|3.ip - (|3.3p has the following property: 



(3.7) P>2ei uniformly in e. 

Namely, this is clear since the function v — p — 2 1\ solves the problem 

v t + (uv) x = ev xx , «|t=o>0. 

From now on, we assume that the initial data p e and Uq belong to the Sobolev space H 2+ @(W.) 
for some < (3 < 1 and satisfy 

< Po < M, HugHoo < Ul, 

and 

ul-nto, Po -> po in Lj oc (M.), 
where ui := |ju ||oo and M := ||po||oo- 

Lemma 3.2. Let (u, p) be a smooth bounded solution of the Cauchy problem (|3.1[) - (|3.3p . Then 
there exist positive constants c\, pi, W± , and Z\ such that 

2e 1 <p<p l , |m| :=p\u\ <c lP (1 + | lnp|) <m l} pi := (2 e% + M) e Ul , 

mi := ci sup p(l + | In p| ) , 
(3 8) °<p<pi 

< W := pe u < Wi, 0<Z := pe~ u < Z u 

uniformly in e. 

Proof. Passing to the Riemann invariant variables 

w := u + In p, z := u — In p, 
we can rewrite the system (|3.1[) - (|3.2|) as 



wt + w x [ u + 1 h — — ) = ew xx — 



2ei _ 




3ew- 


P 




4 


2ei _ 


ew x 


3ez r 


P 


2 


+ 4 



#t + Zx(li — 1 H — : j — • — ' -., : . 

p 2 4 4 

By the maximum principle, 

w < maxrao(i), 2 > minzo(a;). 
Now, the estimates (3.8) is a simple consequence of these inequalities. □ 

By the estimates (3.8) there exist sequences W e ™ , Z e ™, p e ", and to 6 ™ := p en u e ™ and a family of 
non-negative probability measures v x ,ti called Young measures, defined on the (W, Z)-plane, such 
that 

(3.9) W e " W, Z, p e - p, p e -M e " m weakly * in L£ C (II), 

and 

(F(VT"(:E,t),Z e "(x,t)) - (F)) (p{x,t)dtdx -> 0, 

where we have set (F) := J™^ F(W, Z)dv x ^ for any test function 95 £ D(K 2 ) and any continuous 
function F(W,Z) e C ioc (R 2 ). Moreover, 

suppi/ a ,, t c{(W,Z):0< W< Wl, 0<Z<^i}. 

For a proof that to each bounded sequence v n (x,t) one can associate a Young measure /x X) t, we 
refer to Tartar [32] and Ball pQ; see also [30] . 



cS 



LEFLOCH AND SHELUKHIN 



Lemma 3.3. (Entropy dissipation estimate.) The smooth solution (u, p) of the Cauchy problem 
(|3.1[) - (|3.3[) satisfies the estimate 

(3.10) fA + efm l\\ Llc{n) < c 

uniformly in e. 
Proof. The identity 



, , (1 + plnp-p)) + ^ + epu 2 x 

(3.11) dtK 2 7/3 

" fP U , 1 1 1 ^ ClW 1 T 

= ~~q^{~ + u /° ln P - eftrlnp - 2eiulnp - e( — )x ^— } =: -J a 

follows immediately from (|3.ip and (|3.2[) . Multiplying this identity by the function ty(x) intro- 
duced in the proof of Lemma 3.1 and integrating with respect to x we deduce, in view of the 
estimates (3.7) and (3.8), 

J JV x dx<~J + epu 2 x )dx + cj^(l + ^-)dx. 

Hence, we have 

J J *^ +tpul)dxdt < c, 
which yields the desired estimate. □ 
We rewrite the equations (3.1)-(3.2) as a quasi-linear parabolic system: 

(3.12) u t + ai(u, p, u x p x ) = eu XXl p t + a 2 {u, p,u x p x ) = e p xx , 
where we have set 

p x 2ep x u x 2eip x 

oi := uu x 5—, a 2 := (pu) x - 2e\ u x . 

P P P 1 

In view of (3.7) and (3.8), we obtain the global a priori estimates 

2e%<p<pi, \u\ < c(ui,px,ei). 

With these estimates at hand, it is a standard matter to derive estimates in Holder's norms, 
depending on e, by standard techniques of the theory quasi- linear parabolic equations [13j . We 
will only sketch the derivation. Let C( a; j^) be a smooth function such that £ 7^ only if x £ u>, 
where u> is an interval [xq — a, xq + a] . Denote 

ir 71 ' := max{ it — n, 0}. 

Multiplying the second equation in (|3.12p by £ 2 p^ and integrating with respect to x, one obtains 

j t J e\p {n) \ 2 dx + e J C 2 \p x n) \ 2 dx< 1 J(C + C\Ct\)\p {n) \ 2 dx + 1 J (l p > n dx. 

Similarly, for the velocity variable one gets 

| J C 2 \u^\ 2 dx + e J C>W| 2 dx< 7 y - (C| + ClC t |)h ( " ) | 2 + Clp>n+eC 2 |p£" ) | 2 ^. 

These inequalities imply that u and p belong to a class r E>2(Q, M, 7, r, S, n) [13] (Chapter II, §7, 
formula (7.5)), for some parameters Q, M, 7, r, 5, and n. Then it follows that the estimate 

\\u,p\\ H <*.°'/2(ux[o,T]) < c 

holds for some a £ (0, 1). 

In the same manner, one can estimate the Holder norm of the derivatives u x , u XXl u±, p x , p xx , 
and pt, in the same way as done in [TT] for a general class of parabolic systems. 

We now arrive at the main existence result, concerning the viscous approximation (3.1)-(3.3). 
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Lemma 3.4. (Existence of smooth solution of the regularized system.) Let Uq, p% £ L°° n H^ oc , 
< f3 < 1. TTien £/ie Cauchy problem (3.1)-(3.3) has a unique solution such that 

u,peL°°(n)nHl + /' 1+0/2 (n). 

Now, we set e\ = e r , r > 1, and study compactness of the viscous solutions (w e , p £ ) when e — * 0. 
Lemma 3.5. Given an entropy entropy-flux pair (r)(m, p),q(m, p)), m = pu, the sequence 

dt dx 

is compact in WZ C ' (II), where rf = rj(m IL ,p' i ), q e = q(m e ,p e ). 
Proof. We use the following lemma due to Murat's lemma |23j . 

Let QcK 2 be a bounded domain, Q £ C ' . Let A be a compact set in W~ 1,2 (Q), B be a bounded 
set in the space of bounded Radon measures M(Q), and C be a bounded set in W ' P (Q) for some 
p £ (2, oo]. Further, let D C D'(Q) be such that 

D c {A + B) n C. 

Then there exists E , a compact set in W~ ' iff) such that D c E. 
By definition, the functions rj(m, p) and q(m, p) solve the system 

2m to 2 

P P 2 

Hence, calculations show that 

(3.13) F = 2e x m x { V l + ^) + 2 Cl (-^ - ^ + ^) + + erf m m ra = 

p /r p z p a p H 

£i^(<?m + >7p) - 2e 1 ^^ + erf xx - e\rf pp p\ + rf mm to 2 + 2??^ m p a .TO a; ]. 

We check the conditions of Murat's lemma. By Lemma 3.2, the sequence 9 e is bounded in 
W ; ~ c '°°(n). Hence, it is enough to show that er^ — > in L 2 oc (n) and the residual sequence 
° e - vrilx is bounded in L\ JJl). 
We have 

<f m + vp 



Thus, by estimates (3.8) and (3.10), erf x -> in L 2 



&) x = epu x r/ m + ep x - 

2 

loc' 



Consider the sequence 8 e — ef]% x . We have 

° e ~ rfxx = -4vp P pl + V e mm ml + 2rf p x vn x \ + e 1 u x {q t m + if ) - 2e x -^ 



Ipplx 1 Irani x 1 Ipml- 

Each term on the right hand-side is bounded in L\ oc provided t\ = e. Indeed, by (3.7), 

■,V2|„.J _ , IrAnJ -/2l\p x \ 



2e 1 p 1 ' 2 \u x \ r— 1/2 2e x \p x \ 

2ei\u x \ = m < V2ep ' \u x \, < 

P 1 ' 2 P r ■ 

Moreover, if t\ = 0(e), 

(3.14) eiUx{qk + r] e ) _ 2€l P^k^ in L 2 o(n) . 

P 

The other terms are treated similarly. This completes the proof. □ 
Given two entropy pairs {rji(m, p), qi(m, p)), (i = 1,2), from Lemma 3.5, we define 
fk(W,Z)=T h (f 2 (W,Z),f 1 (W,Z)), q i (W,Z) = q i (f 2 (W,Z),f 1 (W,Z)). 
Clearly, the functions 

dtvt + dxqt 

are compact in W,~ ' (n). Hence, by the div-curl lemma [32] . Tartar's commutation relation 
(3-15) (fa q 2 - f] 2 qi) = (%) (q 2 ) ~ (fj 2 ) (q~i) 
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is valid. 

For reader's convenience, we remind that the div-curl lemma states the following. 
Let Qcl 2 be a bounded domain, Q € C ' . Let 

w 1 ^w, W2 -»■ W2, «i — >• vi, «2 v 2 , 

weakly in L 2 {Q), as k — > 00. With curl(w\,W2) denoting dw2/ 'dx\— dw\/dx2, suppose that 
the sequences du^ff, i;^) and curl(wi, lie in a compact subset E of W~ 1,2 {Q). Then, for a 
subsequence, 

v\ w\ + V2 u>2 — > V\W\ + v 2 w 2 otD'(Q) as k — > 00. 

The further claim is due to the fact that system (|1.1|) - (|1.3|) is invariant with respect to the 
scaling p — > Xp. 

Lemma 3.6. If (m, p) is an entropy solution with initial data (mo,po) then (cm,cp) is also the 
entropy solution with the initial data (cmo, cpo) > where c is an arbitrary positive constant. 

Proof. The claim follows easily from the fact that the pair (rj(cm, cp), q(cm, cp)) is an entropy- 
entropy flux pair as soon as the pair (r;(m, p), q{m, p)) is an entropy-entropy flux pair. 
Given A > 0, let us consider the auxiliary problem 

(3.16) p t + {pu) x = ep xx + 2e 2 u x , 

(3.17) (pu) t + (pv?) x + p x = e(pu) xx + e 2 (u 2 ) x + 2e 2 (lnp) x , 

(3.18) p\ t =o = \p e Q {x) + 2e 2 , «|t=o = u e ( x )i 
where €2 = Aei = Ae r . 

The main feature of the auxiliary problem is the following. If the functions (u £ , p e ) solve the 
problem (3.1)-(3.3) then the functions (u e ,p e ) solve the problem (|3.16[) - (|3.18[) with p t = \p e . 
The solution (it e , p e ) of problem (|3.16|) - (|3.18[) obeys the estimates 

(3.19) 2e 2 < p e < (2e 2 + A|| j0o ||oo)e ll " olU =: p 2 , \u e p e \ < cp e (l + | lnp e |) 

uniformly in e. Lemmas 3.3-3.5 are also valid for (u e , p e ). The corresponding Young measure v x t 
has a finite support: 

(3.20) supp v x>t C {(W, Z) : < W < W 2 , < Z < Z 2 } := K. 
We impose the following smallness conditions for A: 

(3.21) p 2 <l, In— >A. 

P2 15 

Assume that the solution (it e , p c ) of the auxiliary problem converges to an entropy solution (m, p) 
of the problem (|2.ip : 

(u e p e , p e ) — > (m,p) almost everywhere in LT. 
The initial data for (m, p) are 

p|t=o = Apo, m| t=0 = Am . 

By Lemma 3.6, the functions (m , p ) = (m/X, p/X) is an entropy solution of the same problem 
with the initial data 

p |t=o = po, |t=o = m . 
Thus, it is enough to study convergence of the solutions to the auxiliary problem. 
With the condition (|3.21[) at hand, the function 

D(R):=(-± + ^B)eW, R:=lnp, 

from Section 5 admits the estimate D(R) > ^e R ^ 2 . Hence, D(R) vanishes only at the vacuum 
points p — 0. 
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To conclude the section, we remark that the parameter t\ serves as a regularizer for the hyper- 
bolic system ([rT j) -(fO ]l with e = due to the estimate ((XT)) (cf. [2"T)). 

4. A LARGE CLASS OF MATHEMATICAL ENTROPIES 

4.1. Symmetry group analysis. We already pointed out that a pair (r],q) is a mathematical 
entropy if and only if rj satisfies 

(4.1) T)pp = -^Vuu- 

In order to derive an explicit formula for the weak entropies of the Euler system we rely on 
symmetry group analysis, following [31]. Using the Riemann invariants 

w := u + In p, z := u — In p, 

we write the equation (14.1)1 for the entropies in the form 



(4.2) F(rj w ,ri z ,rj wz ) := r\ wz - A(r] z - rj w ) = 0, A := i. 

In the more general case where A is a function of w and z, complete group analysis arguments 
were developed in Ovsyannikov's monograph [28] . In our case, A is a constant and this analysis 
is simpler. We only present the results of the formal derivation and we refer to [28) for further 
details on the theory. 

We look for a one-parameter group determined by the infinitesimal operator 

X =£(w,z,r))-^ +t(w 1 z 1 t 1 )^- +ip(w,z,ri)-^-. 

Calculation of the first and the second prolongations of this operator yields 

fS fS Q Q £\ 

X 1 = X + C w — h C* -5 — , x 2 = x 1 +C wa ~ — +C Vu "^ — + C" 



where 



drjio dr\z drj ww dr\ wz dr\ z 

d d 

C w = D w ip - r) w D w £ - t] z D w t, D w = - h j]. 



d d 

C" =D z cp-r) w D z Z-r) z D z T, D z = —+r] z — 



w or] 
d_ d_ 

Oz oi] 



and 



C Vwz =D z ip w + ■qwD^n + VnVwz ~ VwwD z ^ - r) wz (D w £ + D z t) 

- riw{Dz£,w + VwDz&j + £, v Vwz) ~ VzzD w t - t] z (D z t w + r] w D z T n + r) wz T n ). 

Note that we need not calculate the coefficients and C v " ■ Application of the operator X 2 to 
F and analysis of this application on the manifold F = enable us to conclude that the equation 
(|4.2[) admits four one-dimensional groups G; and one infinite-dimensional group G5, associated 
with the infinitesimal operators 

d d d d d d d 

^— , V^~, E: = w 7i z— + A(w + z)r]—, P(w,z) — 7 

ow oz or] ow oz or] or] 

where (3 is a solution to (|4.2j) . The fact that the equation (|4.2p admits the group Gi means the 
following: if i](w, z) solves (|4.2I) then for any c, £ € i? the following functions are solutions of (|4.2[) 
as well: 

r](w + c, z), r](w,z + c), cr](w,z), 
r](e~^w, e^z) exp (Aw(l — e~^) — Az(l — e^)), f]{w, z) + (3(w, z). 

Note that, once this assertion is obtained, its validity can be checked directly without referring to 
group analysis. 

Let us find an invariant solution to the equation (|4.2[) , by applying the one-dimensional group 
G4 associated with the infinitesimal operator E. First, we look for invariants I(w, z, rj) of the 
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group G4 as solutions of the equation E(I) = 0. By the method of characteristics, one derives 
the system of O.D.E.'s 

dw dz drj 
w z A(w + z) 

and obtains easily the following two invariants: 

h=wz, h=r,e~ A[w ' z) . 

Then, we look for an invariant solution of equation (|4.2D in the form (see again [28]) I2 = f{Ii), 
or equivalently 

(4.3) n = e Mw ~ z) f{wz). 

Plugging (|4.3[) in (14. 2\i , we arrive at the following condition for the function f(x): 

(4.4) xf"(x) + f(x)+A 2 f(x)=0. 
In conclusion, the equation (|4. 2|) admits the solution 

r 1 = p 1 ' 2 f(u 2 - In 2 p), 
where the function / satisfies the equation 



4.2. Mathematical entropies. We search for entropies 77 = n(p,u) having the form 

7]{p,u) = p 1/2 f(u 2 -In 2 p). 

It is straightforward to see that 77 solves the entropy equation (14. ip if and only if the function 
/ = f(jn) is a solution of the ODE 

(4.5) mr + f + A 2 f = 0, A 2 = ±- 

16 

With the notation 

R := In p 

the entropy then takes the form 

r ] = V (R,u) = e R / 2 f(u 2 -R 2 ), 
while the entropy equation (|4.ip reads 

(4.6) L(r}) := T] RR - T] uu - iiR = 0. 

One solution to the second-order equation (|4.5|) is given by the following expansion series 



^2, A n o 

f( m ) : =E(^r) (-1)^ 

with 

" f A n y r 



/(0) = 1, f(0) = ~A 2 , fhy 2 ) = Y,( 



n=0 



Observe that /(m) can be represented by the Bessel function of zero order. Given any function 
g : K — > M, we introduce the notation 



g(m) 



g(m), m < 0, 
0, m > 0. 



In particular, this defines the function /. We denote by S the Dirac measure concentrated at the 
point m = and, more generally, by S m=a the Dirac measure concentrated at the point a. We 
denote by 2)(R) the space of smooth functions with compact support and by D'(1R) the space of 
distributions. 

Lemma 4.1. The function f solves the ordinary differential equation (|4.5p in D'(R). 
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Proof. Given a test function (p 6 D(R), we have 

(mj",<p):= [(mcp)"fdm= f (mip)" f dm = (/(0) 6 + ^hj 77 , <p) 

and 

Hence, we find 



(mf +f + A z f, v ) = (m f" + f + A 2 /, ip) = 0. 



Motivated by Lemma 4.1 we introduce the function 

X {R, u) := e«/ 2 J(u 2 - R 2 ) = e R ' 2 l R *_ u *> f(R 2 



□ 



u 2 ) 



(4 ' 7) = ^ 

where 

A + :={ A ' 

[0, A<0, 

and l g >o denotes the characteristic function of the set {g > 0}. 

Theorem 4.2. (Existence of the entropy kernl.) The function \ defined by (|4.7|) is a fundamental 
solution of the equation (14. 6|) in D'(R 2 ). More precisely, we have L(x) = iSr RiU -\—(o t o). 

Proof. From the definition (|4. Tjl of \ and since the multiplicative factor e R / 2 is a smooth function, 
it is straightforward to obtain 

L(x) = e * /2 g RR -7 uu -\) 

in the sense of distributions. Note that, throughout the calculation, / = f(u 2 — R 2 ) and that the 
variables u and R describe M. We compute each term in the right-hand side of the above identity 
successively. We have first 

(fRR><P) = (fiVRR) = // fip RR dudR 

J J u 2 — R 2 <0 

f ip RR dRdu 

\R\>\u\ 



II (<p R f) R + 2Rf ifRdRdu. 

J J\R\>\u\ 



Hence, we obtain 
(7rr, V> 



/(0) / ((pR(-\u\,u)-<pR(\u\,u))du+ // ({2Rtpf') R -tp(2f ~AR 2 f"))dRdu 

l\R\>\u\ 



/(Q) / (MH«U)-MH*0)^-2/'(o) / (v>(-|«|,«)H + ^(H,«)H)du 

<p (4i? 2 /" - 2/0 dfldu. 



|i?|>|u| 



Thus, we have established that 



(4.8) f RR = AR?f»-2f>+J 1 , 
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where J\ is the distribution defined by 

r+oo 



{Ji,V>) =/(0) ((p R (u,u) - tp R (-u,u)) du + J (ip R (-u,u) - (p R (u,u)) duj 

+ 2/'(0)^ / (u(p(u,u) + uip{—u,u)) du — I it, u) + uLp{u,uf) duj . 



The derivative in u is determined in a completely similar fashion. We get 

(/««></>) = / / fifiuududR 

JJ\u\<\R\ 
*+oc /. + |-R| 

/ ((fuf)u - 2u ip u /') dwdi? 

-oo J-|-R| 

+00 

=/(0) / (^„(i?,|i2|)-^(i?,-|«l))^ 



(2/' + 4u 2 f") -2(uf <p) u ) dudR 

\u\<\R\ 



and thus 



(4.9) f uu = 2fi + 4u?f" + J 2 , 

where the distribution Ji is given by 



(J 2 ,<p) = f(0)(J v u {R,-R)-v u {R,R)dR + j ^ u {R,R)-^ u (R,-R)dR^j 

. rO r+00 

+2/ , (0)(/ (Rip{R,-R) + R<p(R,R))dR- (Rf(R, R) + Rip(R,-R)) dR 
^ J-00 Jo 



Now, since the function / satisfies the differential equation (|4.5p it follows from (|4.8[) - f)4.Q() that 

f RR - f un - //4 = Jl - ^2- 

To conclude, we observe that 

/0 />oo 
uip{—u,u)du=— I Rip(R,—R)dR } 
-00 Jo 
/•oo />0 

/ uip(—u,u)du=— / Rip(R,—R)dR, 

Jo J -co 

and 

/o <-o 
ip u {R,R)dR = ip(0) - / ip R (R,R)dR, 
-OO •/ —OO 

/•OO /"OC 

/ <^(.R,i?)di? = -^(0)- / f R (R, R) dR, 
Jo Jo 

/0 <-0 
^(ii,-iJ)dR = -^(0)+ / <p R (R,-R)dR, 
-OO J —OO 

poc poo 

/ -iZ)dH = p(0) + / <p R (R,~R)dR, 

Jo Jo 

we find that 

(Ji-J 2 ,V>)= 4/(0)^(0). 
Since /(0) = 1 and e fl / 2 = 1 when R — 0, this completes the proof of Lemma 4.2. □ 
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Lemma 4.3. The kernel x vanishes on the vacuum 

lim x(i?, u) = for every u, 

R.-^ — ao 

and, at the origin R — 0, satisfies 

(4.10) |m X (R,-) =0, M Xfl (fl ) -)-»±2d B=0 

m £/ie distributional sense in u. Moreover, for any fixed R, x has a compact support, precisely 

X(R,u) = 0, \u\>R. 

It is smooth everywhere except along the boundary of its support where it has a jump of strength 
±e R l 2 . 

Proof. Detailed behavior of x as R — > — oo can be derived from the asymptotic formula [25) 

~ ,x n ^ e 2x , „,l 
,1 

o 



£(-) 2 = VT^^ 1 + °(-)) when x T °°- 

It follows that 



e (-\R\+VK 2 -u->)/2 1 

X (R, u) = l fl 2_ M 2> ^ {R2 _ u2)1/i (1 + 0(-j===)) when R | -oo. 
Next, given <p = (f(u), ip — ip(R) € D(M) we have 

(XR,V^} = - / <p(u) / e R/2 f(m)ip R dRdu 

JR J\R\>\u\ 



J+ I ipil>e R/2 (| -2Rf')dudR, 



where 



l«l>M 

| ti | /"OC 



J= / p(u){ / + / X(e R ' 2 f^) R dRdu 
p(u) (e-l"!/ 2 - e |u|/2 d- 



(p(w)e-l u l/ 2 V(-|M|)^= / e R / 2 ip(R) (<p(R) + <p(-R)) 1 R<0 dR, 



We calculate 



V(u)e^/ 2 i;(\u\)du= f e R / 2 4>(R)(^(R) + ^(-R))lR> dR. 
It follows that, for each R, xr is a distribution in the variable u, given by the formula 
(XR{R,-)M«))= I v {u)e R i 2 ( f -^—^±-2Rf'{u 2 -R 2 ))du 

J\u\<\R\ \ Z > 

+ e R/2 (tp(R) + <p(-R)) (1r>o - 1r<o). 
This completes the proof of (|4.10[) and thus the proof of Lemma 4.3. □ 

Since the equation (|4.6[) is invariant under the transformations u i— ► u — s for every constant s, 
we deduce immediately from Lemma 4.2 that, for every s € K, the function 

x(R,u-s) = e R / 2 7(\u~s\ 2 -R 2 ) 

satisfies the partial differential equation 

(4.11) L( X )(R,u-s) = 4:6 iRjU)={0tS) 

in D'(M 2 ). We arrive at : 
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Theorem 4.4. (The class of weak entropies to the isothermal Euler equations.) Restrict attention 
to the region R < (respectively, R > 0). The formula 

r}(R,u)= / x(R,u- s)ip{s)ds, 
Jm. 

where tp is an arbitrary function in i 1 (R) describes the class of all weak entropies to the Euler 
equations for isothermal fluids (1.1) —(1.3) . In particular, for all u el we have 

(4.12) lim T)(R,u) = 0, lim r, R (R,u) = ±2ip(u), lim r](R,u) = 0. 

R.^0 il— >0± R—y—oo 



Proof. It follows from (|4~TTj) that, for all <p £ D(R 2 ), 

/ L(rf)<p dRdu = 4 / ip(s)ip(s,Q)ds, 
it Js. 

which implies that 

L(r/)=0, R^O. 

Since, for any fixed s,R, the fundamental solution \{Rt u ~ s ) nas a compact support in the 
variable u, we also have 

I x(R,u-s)ip(s)ds-+0, R^O. 
Js. 

□ 

4.3. Mathematical entropy-flux functions. We look for the entropy-flux kernel a which 
should generate the class of entropy flux-functions q via the general formula 



q(R,u) — J a(R,u, s) ip(s) ds. 

In the variables (R,u), the system of equations characterizing the entropies 

q P = un p + p- 1 rj u , q u = pr) p + ur) u 
reads, by setting Q := q — ur), 

(4.13) Q R = r) u , Qu = VR~ il- 
lt is clear that the entropy flux can be deduced from the entropy by integration in R and u. We 
focus attention on the region < p < 1, that is, R < 0. We will use the notation 

a V b := max(a, b). 

Theorem 4.5. (Entropy-flux kernel.) The entropy flux kernel has the form 

<j(R, u,s) = u x(iZ, u — s) + h(R, u — s), 
where the function h admits the following representation formulas: 

d r R 

h = —sgn(u — s) + — — / x( r > u ~~ s ) dr, 
ou J 

or equivalently 

(4.14) h=-^-H(\u-s\,R), H=\u-s\+ f e r / 2 f(\u - s\ 2 - r 2 ) dr, 

" s J -(\R\V\u-s\) 

or still 

h =sgn(u - s) ( e -l"- s l/ 2 l| u _ a |<| fi | - 1) 

(4.15) 

' ' (u-s)e r/2 f'(\u-s\ 2 -r 2 )dr. 
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Proof. In view of (|4. 13|) we can calculate any value = Q(R*,u*) via an integral, as follows 

Q* = J Vu dR + (t]r - rj) du, I* = l~ U Iq, 

where l~ and Iq are the curves in the (R, w)-plane given by 

l~ : R = 0, u = X, — oo < A < u*, 
Iq:R = AR*, u = u*, < A < 1. 

It follows from [g7Ig| that 

ip(u)du+ / r) u (R,u*)dR. 

-oo JO 

Substituting Z~ by 

Z + : u = A, < A < oo, R = 0, 

one obtains similarly that 

Q* = -2 / ijj(u)du+ / rj u (u*,R)dR. 

Joo JO 

Observe, that 

V>( u ) du + i>{u) du = j ijj{u)sgn{u* — u) du. 

-oo Joo Js. 

Hence, 



Q{R,u) — — / tp(s)sga(u — s) ds + j ip( s ) q- / x( r , u — s ) drds, 



Next, we have 

j X(r,u~ s)dr = - / e r/2 f(\u - s\ 2 - r 2 )l r< _ ]u _ s] l r> _ ]R] dr = -Hi, 
Jo J-\R\ 

where Hi is the last integral in (|4.14[) and, therefore, the first formula is established. 
To calculate 

Q CJ /— \U-S\ 

4rZi = Zrl e^ 2 f(\u-s\ 2 ~r 2 )dr, 

l-(\R\V\u-s\) 



we observe that 



du du 
8 



du 



(\R\ V \u - s\) = l| u _ s |>|ij| sgn(u - s). 



Hence, we have 

■^-Hi=2 [ ' (u-s)e r/2 f'{\u-s\ 2 -r 2 )dr-f(0)e-\ u - s V 2 sgn(u-s) 

+ l\u-s\>\R\e- mvlu - sl)/2 f(\u - s\ 2 - (\R\ V \u - s\) 2 ) sgn( U - s). 
The last term coincides with 

l\u-s\>\R\ e- |u - s|/2 sgn(u-s). 
Thus, the representation formula (|4. 1 5[) is proved and the proof of Theorem 4.5 is completed. □ 
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4.4. Singularities of entropy and entropy-flux kernels. From the above results we see that 
the functions x an d h are continuous everywhere except along the boundary of their support, 
that is, the lines u = s ± \R\. The most singular parts (measures and BV part) of the first order 
derivatives of the functions x an d h with respect to the variable s are now computed. 

Theorem 4.6. (Singularities of the entropy kernels.) The derivatives Xs and h s in D'(R) are as 
follows: 

(4-16) Xs = e R/2 (Ss=v,-\R\ ~ <^=«+|fl|) + G x {R,u~ s) l|^_ s | < | jR | , 



(4.17) h s =e R/2 [8 s=u _\ m +5 s=u+lRl j + G h {R, u - s) l|„_ s |<|fl|, 

where, for all \v\ < \R\, 

G x (R,v) :=2e R/2 vf'(v 2 - R 2 ) 

and 

G h (R, v) := e-l-l/ 2 (1/2 - 2\v\) + 2 / (e r ' 2 f'(v 2 - r 2 ) + 2e r ' 2 v 2 f"(v 2 - r 2 )) 

J-(\R\V\v\) V ' 



dr. 



It will be convenient to extend the functions G x and G h by continuity outside the region 
\v\ < \R\ by setting 



G x {R,v) 

and 



2\R\f'(0)e R / 2 , v>\R\, 
-2\R\f{Q)e R ' 2 , v<-\R\, 



G h {R,v) = e R ' 2 (2R+l/2) 1 \v\>\R\. 
Proof. Given a test function ip = ip(s), we can write 

ru+\R\ 

X tp'(s)ds = e R / 2 / ip'{s)f(\u-s\ 2 ~R 2 )ds 

Ju-\R\ 

u+\R\ 

= e R / 2 (<p(u+\R\)~cp(u- \R\)) +e R ' 2 J 2<pf (\u - s\ 2 - R 2 ){u - s)ds, 

u-\R\ 

which yields the first formula (|4.16|) . 
Next, it follows from fl~T4"]) that 

h s = sgn(u - s) fH"-^/ 2 (s s=u _\ R \ - 5 s=u+ m + ^sgn(w - s)l|„_ s |<|ij| 

-|ll-s| Q 

^- ( (u - s) e r/2 f'(\u - s\ 2 -r 2 ))dr + 2 /'(0) e -^ /2 (u - s) sgn(s - u) 



-2 



-2 e-(l*l v l»-D/ a l {u - sl > lR{ f(\u - s\ 2 - (\R\ V\u- s\) 2 ) (u - s) sgn( W - s). 
The last term above coincides with 

-2e-l"- s l/ 2 l| a _ s |>| fl | /'(0)|u-*| 
and, therefore, the second formula (|4.17[) is also established. □ 
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5. Reduction of the support of the Young measure 

5.1. Tartar's commutation relations. We now turn to investigating Tartar's commutation 
relation for Young measures, following the approach in Chen and LeFloch [4J [5]. In the previous 
section we constructed the class of weak entropies r\ and entropy fluxes q in terms of the variables 
p and u. We can also express rj and q as functions of the Riemann invariants W and Z , via the 
following change of variables 

fj(W,Z);=r)(u,p), q(W,Z)=q(u, P ), 
W := pe u , Z := pe~ u . 

To simplify notations, it is convenient to adopt the following convention. In the rest of this section 
we will write (F) — J F(u, p) dv instead of J F(W, Z) dv, by assuming that p, u are the functions 
of the variables W, Z given by 

p={WZf'\ u=\\n^. 

We will prove: 

Theorem 5.1. (Reduction of the support of the Young measure.) Letv = v{W,Z) be a probability 
measure with support included in the region 

{{W,Z) : < W < W 2 , 0<Z<Z 2 } 

and such that 

(5-1) (77a q 2 - T] 2 qi) = (in) (q 2 ) - (m) (qi) 

(where (F) := (v, F) ) for any two weak entropy pairs (rji, qi) and (772, q 2 ) of the Euler equations 
(1.1)-(1.2). Then, the support of v in the (W, Z)-plane is either a single point or a subset of the 
vacuum line {p = 0} = {WZ = 0}. 

The proof of Theorem 5.1 will be based on cancellation properties associated with the entropy 
and entropy-flux pairs of systems of conservation laws. The key idea (going back to DiPerna [10] ) 
is that, nearby the diagonal {s 2 = S3}, the function 

E(p, v; s 2 , s 3 ) := \(p, v - s 2 ) <r(p, v, s 3 ) - x(p, v - s 2 ) cr(p, v, s 3 ) 

is much more regular than the kernels x an d o~ themselves. 

The principal scheme can be explained as follows. Given functions ipi € D(R), (i = 1, 2, 3), we 
define the entropy pairs 

r]i(u,R) = J x(u- Si,R)tpi(si)dsi, q l {u 1 R) = J a(u,R, 

and deduce from Tartar's relations (|5.ip the following remarkable identity (see Chen and LeFloch 
[I], as well as the earlier work [H)] ) 

(V1Q2 - qim) (m) + ((HV3 - V1Q3) (%) + (q3m - 773*72) (vi) = 0. 

Next, substituting ipi(si) with —tp'^Si) = —Piipi(si) and denoting Fi = F(u, R, Sj), we arrive, 
after cancellation of the arbitrary functions ipi(si), at the equality 

^ ^ (Xi p 2h 2 - h 1 P 2 x 2 ){P3X3) + (hi p 3X3-XiF3h 3 )(P 2 x 2 ) 

which is valid in D'(R) 3 . In view of the expression of the distributional derivative of a and h 
(Theorem 4.6), each term in (I5.2[) can be calculated explicitly. Denoting 

w = u — \R\, z = u + \R\, 

we find 

(5.3) xi^/7.2 - h\P 2 x 2 = e R / 2 {h x - xi)bs 2=w - e m (h x + X i)S S2 = z + (hiG% - xiG£)l| u _. a |<| fl | 
and, similarly, 

(5.4) Xl P 3 h 3 - h x P 3X 3 = e fl/2 (/ii - Xi)S S3=w - e*' 2 ^ + Xl )6 S3 = z + faGg - Xi^)1| u _ S3 |<|h|. 
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Moreover, we have 



- F 3 /i 3 P 2 X2 =2e R (8 S2=z 5s 3=w - S S2 =w S S3=Z 



+ e R / 2 (s s3=w (G 2 - G*) - 5 S3=Z (G% + G%yjl\ v -sa\<\R\ 



+ (G3G2 - G|G3)l| u _ S2 | < | fl |l| M _ S3 | < |i ? |. 

In view of the formulas (|5.3p and (|5.4[) the right-hand side of (|5.2p contains products of functions 
with bounded variation (involving a and h) and Dirac masses plus smoother terms. Such products 
were earlier discussed by Dal Maso, LeFloch, and Murat [7]. On the other hand, the right-hand 
side of (|5.2p is more singular and involves products of measures, product of BV functions by 
measures, and smoother contributions; see (|5.2p . Our calculation below will show that the left- 
hand side of (|5.2p tends to zero in the sense of distributions if s 2 — > si and S3 — » s±, while the 
right-hand side tends to a (possibly) non-trivial limit. 

We test the equality (|5.2p with the function 



(5.5) ip(s) ip 2 (s - s 2 ) (fil(s - s 3 ) := ip(s) ^ ip 2 { S f 3 C ^ ) 

of the variables s = sx, S2, S3, where ip € I>(R) and ipj : R — > M is a molifier such that 
Vi(sj) > 0, / ^j(sj) dsj = 1, supp ^(sj) C (-1, 1). 



5.2. Nonconservative products. To provide testing of equality (|5.2|) by the function (I5.5[) . we 
will need the following technical observations. 

Lemma 5.2. Let %p, F : R — > R otic? / : [a', 6'] — > M be continuous functions. Then, for every 
interval [a, b] C R, tte integral 



I e (a,b,a ,b) :=/ V( s ^(si ~ a ) / F(s 3 ) ^ (si - s 3 ) ds 3 dsi 
ftas £/ie following limit when e — > 

^(o) F(a) /(a) l a ,< a<fc , + B 2 " 3 /(a'+) l a=B , + G 2 ~ 3 /(&'-) w) , 
where A^ 3 '■— B2 3 + C 23 and the coefficients B~ and C~ depend only on the mollifying functions: 

B 2,3-= / V>2(yi) ^3(2/3) dy 3 dyi, C 23 := / <p 2 (yi) ^3(2/3) dy 3 dyy. 

JO J — 00 J — 00 «/ — 00 



Formally the integral I e has the form 

rb rb 

I(a,b,a',b'):= $( 8l ) /(si) 6 Sl=a / F( S3 )£ Sl=S3 . 

^ si— a 7 »/ S3— a 

Lemma 5.2 shows that this term can not be defined in a classical manner and that, by regulariza- 
tion of the Dirac masses, different limits may be obtained, depending the choice of the mollifying 
functions. 

Similarly we have 

Lemma 5.3. Let ip, F : R — > R and f : [a', b'] — * R 6e continuous functions. Then, for every 
interval [a, b] C R, the integral 

r-b r-b 

J c (a,6,a',6') := / ip(si) f(sx) (pfat - b) F(s 3 ) tpl(sx - s 3 ) ds 3 dsi 



/las i/ie following limit when e — > 0: 

V(6) F(6) (A+3 /(6) l Q , <b<6 , + fl+3 /(o'+) l a =a' + C+, /(&'-) l 0= y) , 
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where A^ 3 := B% 3 + 3 and 

poo poo />0 />oo 

B t3 : = / (2/1)^3(2/3)^3^1,^3 := / / ^2(2/1) (p 3 (y3)dy 3 dy 1 . 

Along the same lines we have also: 

Lemma 5.4. Let ^/>,-F : K — > K fee continuous functions and let the function f : R — > R &e 
continuous everywhere except possibly at two points a and b with a < b. Then, for every real a 
the integral 

K e (a,b,a):= / V(*i) /(s x ) ip%(s x - a) I F(s 2 ) ip 2 (s x - s 2 ) ds 2 dsi 

JR J a 

has the following limit when e — > 0: 
1>{a) F(a) (/(a) l a<a<b + (C^ /(a-) + B£ 3 /(a+)) l a=a + (C+ 3 /(&-) + B+ 3 /(&+)) l a=6 ). 
We only give the proof of this last statement. Lemma 5.2 and 5.3 can be checked similarly. 

Proof. Making first the change of variables s 2 = s% — ey 2 and then s% = ey± + a, one can write 

r ryi + (a-b)/e 

K e = - / ip{eyi + a) f(eyi + a) (£3(2/1) / F(e(yi - y 2 ) + a) (£2(2/2) dy 2 dy%. 

JR Jyi+{a—a)/e 

Clearly, we have K e — > when a < a or a > 6. 
Now, if a = a we can write 

3 -a ,>fr ~oo />— — - 

K' = Y. K t = -( +/+/ )V>0>i)/(siM(*i - «) /_'„ F( Sl -ey 2 )My2)dy 2 ds 1 . 

^ J —00 Ja Jb J 1 

Consider the first term iff: 



/0 /-3/1 — 

^(ej/i+a)/(ej/i + a)yj3(yi) / F(e(yi - y 2 ) + a)ip 2 (y 2 ) dy 2 dyi, 

-00 J yi 

which satisfies 

(£3(2/1) / ^Pi{y%)dy 2 dy\. 

-00 J — oo 

Similarly, one see that 

r°° ryi 

K^iP(a)F(a)f(a+) (£3(2/1)/ ^2(2/2)^2^1, K\ -> 0. 

JO J-oc 

The other values of a can be studied by the same arguments and this completes the proof of 
Lemma 5.4. □ 

5.3. Proof of Theorem 5.1. 

Step 1. First, we consider the right-hand side of (|5.2p . Let us denote 

dv:=dv(W,Z), du' := dv{W , Z'), w' = u'-\R'\, z' = u' + \R'\. 

Applying the distribution (-P3/2.3-P2X2 — -P3X3-P2/12) (xi) to the test function (|5.5[) . we write the 
integral 



(5.6) / (P 3 /13 (53)^2X2(52) - P3X3(s 3 )P 2 h 2 (s 2 )) (xi{s))4>(s) <Pz(s ~ s 2 ) <£ 3 (s - s 3 ) dsds 2 ds 3 

Jm. 3 

as the sum If, in which, in view of Theorem 4.6, we can distinguish between products of Dirac 



measures 



iP{ X i)(2e R ip%{s - z))<pl(a - «;)) - 2e\%{s - w)) V \{s - z)))ds, 
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products of Dirac measure by functions with bounded variation 

Il~ [ ^(Xi)(e H/2 rt(s-w)) f (G* - G h 3 )<pl(s - s 3 )l lu _ S3l<lRl d S3 )ds 



Hxi){e R/2 ^{s - w)) j (Gj - G 2 > 2 (s - S2 )l| u _ S2 |<| R |d S2 )d S 
— T e — T e 

II ■■= [ ^(xi)(e R/2 rt(s ~ z)) J (Gg + G' 3 l M(s - « 3 )l|u-.,|<|fl| da a ) ds 

^(Xi)(e R,2 vl{s ~ z)) j (G* + G$)<pl(s - * a )l| u _„|<| Jl | ds 2 )ds, 



and a smoother remainder 



!t = / ^(xi)((G^G 2 - G% Gj)l| u _ S3 |<| R | l|„_ S2 | < | i j|^ 2 (s-s 2 )v3(s-S3)dsds 2 rfs3. 



By change of variable we see that the integral 

I l = ~ I I I e R ^(ey + z) X (R\u' - (ey + z)) 



e 

W,Z W'Z 



(va(l/) V3(» + 2|i?|/e) -<P2(V + 2\R\/e) <p 3 {y)) dydvdv' . 
tends to zero : Jf — ► 0. The same is true for the smoothest term J|, in view of the identity 

It = J J j (f G*{R,u,s 3 )ipl{s-s 3 )ds 3 J* G h {R,u,s 2 )<p 2 {s-s 2 )ds 2 

-j G x (R,u,s 2 )<f 2 (s - s 2 )ds 2 j G h (R,u,s 3 )ipl(s - s 3 )ds 3 ) 

ipx(R' i u ' — s ) dsdvdv , 

which clearly tends to 

/ / f^x(Rf,u' -s)l lu _ s{<Rl (GX(R,u,s)G h (R,u, S ) 

- G X (R, u, s) G h (R, u, sfj dsdvdv' = 0. 

We denote 

Q ± :=G x ±G h , Fi = F(R,u,Si), F- — F(R' , u, Sj). 
Let us now consider the term I 2 = I 2 x — I 2 2 in ()5.6j) . We have 

I W = J J J ^ eR/2 Xi ^(s - w) J Q3 ip%(s - s 3 ) ds 3 dsdudv' . 
Therefore, in view of Lemma 5.2, we obtain that I 2 1 tends to 

z R l 2 4>{w)Q-{w)(x\w)\ wl<w<z ,A 2 3 + X \w'+)\ w = w ,B 2 3 



and I 2 2 tends to 



e R / 2 tp{w) Q (w) (x?Mlw'<w<z>A 3 2 + x , (w , +)l w=wl B 3 2 + x ' V ' -)l w=z ,C 3 ^dvdv' ', 
as e — ► 0. We conclude that the limit of I 2 is equal to 

e R/2 V>H Q-(w) (x'H Uo>< w <z> A~ + x'(w'+) l w = w , Br + x'(z'-) l w=z , G") <W, 



where 

A =^2,3^^3,2' =^2.3^^3,2' ^ = C2,3 _( ^3,2- 
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By Lemma 5.3 we can determine similarly that lim e ^ ^3 is equal to 
J I e R ' 2 ,l>{z) Q+(z) (x\z) 1 W > <Z<Z ,A+ + x!{w'+) l z=w , B+ + X '(z'-) l z=z > C + )dvdv' ', 
where 

A + := A 23 — A 32 , B + := B 23 — B 32 , C + :— C 23 — C 32 . 
In conclusion we have identified the limit of the term (|5.10p , it is equal to 

e R / 2 i>(w)Q-(w)(x'(w)lw' <«,<*> A- +x'{w / +)l w=w ,B- +^{z , -)l w=M .C-^dudu f 
R / 2 i>(z)Q+(z)(x / (z)l w , <z<z ,A+ + X '(w'+)1 Z=W ,B+ + X '(z'-)l z = z ,C + )dvdv'. 



+ 



Step 2. We now proceed by studying the two terms in the left-hand side of l|5.2p . We apply the 
distribution (xx{sx)P 2 h 2 (s 2 ) — hi(si)P 2 x 2 (s 2 )} (^3X3( S 3)) to the test function (|5.5p . We write 
the integral 



(Xi (si)P 2 h 2 (s 2 ) - h 1 (s 1 )P 2 X2(s 2 )} (PsXa(s3)) tp(si) <p 2 (si - s 2 ) (p 3 (si - s 3 ) ds 1 ds 2 ds 3 
as the sum 1 Jf , where 

J{ := / (e R/2 (hi - xi) S S2=W ) (P3X3) V>( s ) fl(s - s 2 ) (p 3 (s - s 3 ) dsds 2 ds 3 



J% {e R/2 (hi + Xi)S S2=z ) (P3X3) V-(s) ip%{a - s 2 ) tp\{a - s 3 ) dsds 2 ds 3 , 

JR 3 

and 

^3 ; = / (( h iG 2 - XiG2)l\u-s 2 \<\R\)(P3X3H(s)<p e 2 (s - s 2 )ip e 3 (s - s 3 ) dsds 2 ds 3 . 



The application of the distribution (xi-f^l — ^l-PaXl) (P2X2) to the test function (|5.5|) can be 
represented similarly. The integral 



{Xi(s)Pah 3 (s 3 ) - hi(s)P 3 x 3 (s 3 )} (P 2 X2(s 2 )) ip(s) y%(s - s 2 ) ip%(s - s 3 )dsds 2 ds 3 
is equal to J^i Kf> where 

K i : = / (e R/2 (hi -x\) S S3=W ){P 2 X2) ^p(s) (p 2 (s - s 2 ) ip e 3 (s - s 3 ) dsds 2 ds 3 , 

JR 3 

K e 2 :=- (e R ' 2 (h x + xi) S S3=Z ) (P2X2) ip(s)<p e 2 (s - s 2 ) <f 3 (s - s 3 ) dsds 2 ds 3 , 
Jr 3 

and 

K 3 : = / (Oi^f -Xi^)l| u _ 83 |<| fl |) (PiX2) ifr(s)<p\(s- s 2 )y%(s -s 3 ) dsds 2 ds 3 . 



Some further decomposition of these integral terms will be necessary: 

1 1 

where 

J i,i '■= I (e R/2 (hi - xi)$s 2 =w) (e R/2 S S:i=w ) ip(s!)ip 2 (s - s 2 )ip e 3 (s - s 3 ) dsds 2 ds 3 , 



J{ 2 :=- / (e R/2 (hi - X i)5 S2 = w )(e R/2 8 S3 = z ) ^(s^p^s - s 2 )<f e 3 (s - s 3 ) dsds 2 ds 3 , 
Jr 3 

j i,3 : = / (e fl/2 (/n -Xi)^2=w)(G3 1 |«-s3|<|fl|) ip(s)<p 2 (s - s 2 )ip 3 (s - s 3 ) dsds 2 ds 3 , 
Jr 3 

J 2 a~- I (e R/2 (hi +xi)S S2=z )(e R/2 S S3=w ) ip(s)ip e 2 (s - s 2 )(p e 3 (s - s 3 ) dsds 2 ds 3 , 
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J2.2 ■= / {e R/2 {hi +Xi)S S2 =z)(e R/2 5 S3=z ) 4>(s)p 2 {s - s 2 )tp e 3 (s - s 3 ) dsds 2 ds 3 , 

Jl,3 = - / (e R/2 (hi +Xi)ds 2 =z)(G 3 < l\ u - S3 \ < \ R \) *p(s)ip 2 {s - s 2 )<fil(s - s 3 ) dsds 2 ds 3 , 

J 3.i : = / (0i G 2 -XiG 2 )1\u-s 2 \<\r\){z R/2 5s 3 =w) i){s)w e 2 {s- s 2 )ipl(s- s 3 )dsds 2 ds 3 , 
Jm. 3 

JI2 : =- / ((^1^ - XiG^l^.^Kifl^e^ 2 ^^) ■0(s)¥'2( s - s 2)^3(s-S3)dsds 2 ds 3 , 

J 3,3 : = / -XiG2)l| tt - S2 |<| J R|)(G , 3 1 l«- ;S 3|<| J R|> V>( S M(« -«2)^3( s ~ s 3) dsds 2 ds 3 . 

The terms iff l7 K{ 2 , etc. are defined in a completely analogous fashion. 
We can put Jf j in the form 

J tl =^JJJiP e (M)/2 (to Xi) ^2(^ )M Ll ^ 1 ) dvdv'ds 

=\fff Hs)e {R+R ' )/2 (h - x)(s) Mvi)\>=evi+rvMvi + *^p-)dvdv'd yi . 
A similar representation formula is valid for K\ x . In consequence we find 

JU - K U =\fff ^ Z (H+R,) ' 2 ( h X)(s)\s=ey 1+W 

<Pa(yi + W ™ - ) - ^3(2/1) ^2(2/1 + ^ ^ i^dvdv'dyx. 

Clearly, we have 

^.i " #f,i - 0. 

The terms z and i\| l contain the product of measures or the product of BV-functions and can 
be treated in the same manner. In turn, one obtains 

J{ 2 — K\ 2 — ► 0, J 21 — K 21 — > 0, J| 2 ^ ^2,2 ~~ * 0' ^3,3 — ^3,3 ~~ * 0- 

Let us consider the terms J| l and if^ ; , containing the product of a measure and a BV- function. 
By Lemma 5.4, the term 

J ll3 = JJJ^(s)e R / 2 (h 1 -xi)<pUs-w) G 3 ' (p e 3 (s - s 3 ) ds 3 dsdvdv' 
converges toward 

e R ' 2 ^(w) (h- X ){w) G*\w)(l w' <w<z' + 1 w—w' 

(^2,3 + ^2,3) + lt«=z' (G 2 ^ 3 + B 2 3 )^dvdv' , 

hence, 

lim( Jf i3 - ^ li3 ) 

By the same argument we find that the term 

J 3 ,i = JJJ e* 72 ^* - «/) y (/nG£ - XiG 2 l ) <^ 2 (s - s 2 ) d S2 d s dvdv' 



tends toward 
where 



J J e R '/ 2 ii(w')G x (w')Ldvdi;' - J J e R 'I 2 i>{w') G h {w') S dvdv' , 

h{w') + 1^=™ (/i(w-) C 3 - 2 + h(w+) B- 2 ^j + l w , =x (h(z-) Ct 2 + h(z+) B+ 2 ) 
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and 

S := 1 W<W ' <Z xO') + l w >= w (x(w-) C3 2 + x(w+) £^ 2 ) + l w ,= z \x{z-)C^2 + x(z+) B 3,2j ■ 
Hence, 

■*M " ^3,1 ^ - J J e R ' /2 ^{w') {1 W ,= W M W + l w , =z M z ) dvdv 1 , 

where 

M w := G x {w)(h{w-)C- +h(w+)B-) - G h {w)(x{w-)C~ + X (tfl+)B _ ), 

M z = G x (z) (h{z-)C+ + h(z+)B+\ - G h (z) (x(z-)C+ + X {z+)B+^ . 

One can see that the integrals, containing the functions l w '— w and l w >— z cancel each other. 
In a similar way, we can treat the other terms and arrive at the final equality 

e R ' 2 ip(w) Q- (w) x» A- l w , <w<z , + e R ' 2 *P{z) Q~{z) X '(z) A+ l w , <z<z , dvdv 1 = 0, 

resulting from (|5.2D . 

Step 3. Observing that A + = —A~ and 

/(0) 



Q-(w) = -Q+(z) = e fl / 2 ( - ^ + 2|i?| + 2| J R|/'(0)) =: D(R), 
(5.7) If Y e R ' 2 D(R) (i>(w) X » l w , <w<z , + ip(z) x'(z) l w '<z<z') dvdv' = 0, 



we can write 



where 

+00 s 2 



Y= J J ^2(52)^3(53) - <P2(s3)<P3,(s3)ds 2 ds 3 . 

— 00 —00 

As observed in [5], the functions cf2 and cf3 can easily be chosen in such a way that 7^0. 
In accordance with our notation, the equality (|5.7p means precisely 



J I D(~ \n(WZj) (W Z) 1 / 4 (W Z') 1/4 (/((- In ^ ln(Z'W)) ^(ln W) l w<w< y z , 

+ /(- ln(W'Z) In |r) </>(- InZ) l^<i/ z <i/ z /) <W = 0. 

The conditions (|3.21[) guarantee that |D(i?)| > e R ' 2 /2. Since f(—x 2 ) > 1 and the function ip is 
arbitrary, we conclude from the last equality that 

(WZ) 1/2 f f (W'Z') 1/4 dv(W', Z')dv(W, Z)=0 

w,z J J {w<w}n{z'<i/w} 

and 

/ / (WZ) 1/2 f [ (W'Z') 1/4 dv(W', Z') dv(W, Z) = 0. 

J Jw,z J J{w'<i/z}n{z><z} 

We arrive at the following important claim: whenever 

W*Z*^0, (W*,Z*) e supp v, 

we have 

// {WZf/ 2 dv{W,Z)^Q 
J Jw,z 

and therefore 

// (W'Z') 1/i dv(W',Z')=0, 
(5.8) ' 11 1 



II (w'zr /4 dv(w',z>) = o. 

J J iW'<l/Z*rn\Z'<Z*j 
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We will now conclude from (|5.8[) that the Young measure is a Dirac mass or a measure con- 
centrated at the vacuum. At this stage, it is useful to draw a picture on the W, Z— plane, with 
the W— axis being horizontal. One should draw two hyperbolas WZ = 1 and WZ — p 2 , keeping 
in mind that suppf lies below the hyperbola WZ — p\, where the constant P2 < 1 is defined in 
section 3. The hyperbola WZ = 1 helps to picture the set 

M* := ({0 < w' < w*} n {0 < Z' < 1/W*}) U ({0 < w' < i/z*} n {0 < Z' < z*}), 

a union of two rectangulars. The relations (|5.8p imply that M* does not intersect the support of 
p : 

(5.9) W*Z* y£ and (W*,Z*) S supp v =>■ M* n supp v = 0. 

By construction, the hyperbola VFZ = 1 does not intersect supp v. The inclusion supp v C {p = 0} 
holds if no hyperbola WZ = S, < 6 < 1, intersects supp;/. If supp^ contains a point (W, Z) 
such that p(W, Z) ^ 0, there is a number < S < 1 such that the hyperbola WZ = S, intersects 
supp v. Let < (5o < 1 be the largest number such that the hyperbola W Z = 5q intersects supp v. 
By (|5.9p , the intersection 

supp^ n {VKZ = S } 
may contain only one point (W* , Z*) and 

supp^ n {0 < WZ < 5 a } = 0. 

Thus 

(5.10) v = ch5* + /i, 

with supp /i C {p = 0}. Throughout the paper we use only weak entropies. Hence, putting (|5.10[) 
into the Tartar's commutation relation (|5.ip . we obtain that any two entropy pairs satisfy the 
equality 

(5.11) a(q 2 i]i - qif] 2 ) = a 2 (q 2 r]i - qir) 2 ) 

at the point (W* , Z*). Let us choose the following entropy pair (as in 3 lj ) 

r n = p B *e A ' u , qi = --^-pBi-l e AiU j A . = ^ B .( B . _ 1)j Bl ^ B 2 . 
t>i — 1 

Now, the equality (|5. lip is rewritten as 



a(l - a)p B 1+ B 2 -X e ( Al+ A 2 )u, (^L_ 



Bi / B 



B 2 -l 

Hence, a = or a = 1. This completes the proof of Theorem 5.1. 

6. Convergence and compactness of solutions 

Due to the decomposition (|5. 10|) of the Young measures, the convergence formulas (|3.9p imply 
that 

W e ^W, Z e -± Z, F(W e ,Z e ) F(W,Z) weakly* in L°°(IL), 

for any function F(a,f3), F £ C(K), such that F = at the vacuum set a(3 = 0. (See formula 
(|3.20p for the definition of the compact set K.) Hence, for almost all (x, 

p e := (VTZ 6 ) 1 / 2 -» p = {WZ) 1 ' 2 =: fr{W, Z), 
m £ := (W e Z") 1 ' 2 ln(^) 1 / 2 -» m = (WZ) 1 / 2 ln(^) 1 / 2 =: / 2 (W, Z), 

;= Z e ) l/ 2(ln( ^ ) l/2 ) 2 _ ^ = /3( ^ z) := ^1/2^(^)1/2)3. 

p £ Z f - p z 

Moreover, 

(6.1) F(m e 1 p e ) -> F(m,p) for almost all (i,f)6ll 
for any function F(m,p) such that 

(6.2) F(a,0):=F(f 2 (a,P),Ma,f3))£C(K), F\ afi=0 = 0. 
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Indeed, one can derive the convergence (|6.ip from the following fact: 

v e — > v and (v e ) 2 — > t> 2 weakly in L 2 oc (II) v £ — > u strongly in L 2 oc (n). 

Let us show that (m, p) is an entropy solution of problem (12. ip . To this end we let e and ei go 
to zero in (|3.13j) . (More exactly we should do it in the similar equality relevant to the auxiliary 
approximation.) If functions 77(777, p), q(m, p) obey the restrictions (|6.2p . one obtains 



(77(777%//) - ij(mQ, p e Q ))ip t + q(m e , p e )ip x dxdt -> / (77(771, p) - r/(m a , p ))ip t + q(m, p)ip x dxdt, 

77 (m e , p e )(p xx dxdt — > 



for any <p G D(K 2 ). 

From now on we assume that ei = e r , r > 1. If a function 77(771, p) meets the conditions of 
Theorem 2.1, the derivatives f] m (m, p) and m p (m,p) are continuous on any closed set 

{0<p<pi, |m| <cip(l + |lnp|)}, pi>0, ci>0. 

Hence, by estimate (j3,19[) . 

lei^x^m + ?7p)| = |2eiu :E (^77 m + ?7 P )| < cei(\uu x \ + \u x \) < e 1/2 p 1/2 \u x \(e~ + |u|p 7 e' 5 ), 

where 2j < ,25 = r(l — 27) — 1. Besides, eip _1 |77 m pa;| < ce 1 / 2 ^ -1 / 2 !^^^ 11 ^ 1 . Now, it follows 
from Lemma 3.3 and estimates (|3.19l) that 

eiu x (q m + rip) - 2e 1 r tm p~ 1 p x ^0 in £ 2 oc (II). 

Taking into account the convexity of the function 77(777, p), we send e to zero in (|3. 13|) to deduce 
that the pair (m,p) is an entropy solution of J2.1[) - (|2.2p . The proof of Theorems 2.1 to 2.3 is 
completed. 

We conclude by giving a proof of Theorem 2.4. Let (m n ,p n ) be a sequence of bounded in L°° (II) 
entropy solutions of the problem (|2.ip obeying the restriction of Theorem 2.4. We introduce the 
sequences 

vv n . — pn" j — Pn e 

Clearly, we have 

W n -± W, Z n -»■ Z weakly ★ in L% C (I1), 
and there exist Young measures v x _ t such that, for all F(a,{3) G C; oc (M 2 ), 

F(W n (x,t),Z n {x,t)) -± (v x , t ,F). 
Given two entropy pairs [T)i(m, p),qi(m, p)) from Theorem 2.1, the sequences of measures 
0™ := d t r)i(m n ,Pn) + d x qi(m n , p n ) = d t fj(W n , Z n ) + d x q(W n , Z n ), 

satisfy the conditions of Murat's lemma and are compact in W;~^' 2 (II). (We remind the notation 
q(W, Z) := q(f2(W, Z), fi(W, Z)).) By the div — curl lemma, the Tartar commutation relations 
(|3.15p is valid for the Young measures v x t . Then we argue like in the proof of Theorem 2.1 to 
arrive at the decomposition (15.10[) for v x t . Hence, 

F(W n (x, t), Z n (x, t)) — * F(W(x, t), Z(x, t)) almost everywhere in n 
for any F(a, (3) G C ioc (R 2 ) such that F(a, /3) = Q if a/3 = 0. Denoting 

p=h(W,Z), m = h{W,Z), p n = f 1 (W n ,Z n ), m n = f 2 (W n ,Z n ), 
we pass to the limit, as n — * 00, in the inequality 

r)[m n , p n ) d t tp + q(m n , p n ) d x cp dxdt + J r](m , p ) ip(x,0) dx > 

and check that (to, p) is an entropy solution of the problem (|2.ip - (l2.2p . The proof of Theorem 2.4 
is completed. 
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